Interlayer spin-singlet pairing induced by magnetic interactions 
in an antiferromagnetic superconductor 
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It is shown that interlayer spin-singlet Cooper pairing is induced by magnetic interactions in a 
metallic antiferromagnet of stacked conductive layers in which each layer is ferromagnetically polar- 
ized and they order antiferromagnetically in stacking direction. As a result, the antiferromagnetic 
long-range order and superconductivity coexist at low temperatures. It is shown that Taf > T c 
except for in a very limited parameter region unless Taf = 0, where Taf and T c denote the antifer- 
romagnetic and superconducting transition temperatures, respectively. It is found that the exchange 
field caused by the spontaneous staggered magnetization does not affect superconductivity at all, 
even if it is very large. The resultant superconducting order parameter has a horizontal line node, 
and is isotropic in spin space in spite of the anisotropy of the background magnetic order. We 
discuss the possible relevance of the present mechanism to the antiferromagnetic heavy fermion 
superconductors UPd2Al3 and CePtaSi. 

PACS numbers: 74.20.-z, 74.20.Mn, 74.20.Rp 74.25.Ha 



In this paper, we show that interlayer spin-singlet 
Cooper pairing is induced by magnetic interactions in a 
certain kind of metallic antiferromagnet. We consider a 
layered system of itinerant electrons in which each layer is 
ferromagnetically polarized but the majority-spin alter- 
nates in stacking direction. Therefore, the magetic order 
is characterized by the wave vector Q = (0,0, 7t/c), where 
we have assumed the a and b crystal axes to be parallel to 
the layers, and the c-axis in the stacking direction, and c 
denotes the c-axis lattice constant. It is also shown that 
the exchange field caused by spontaneous staggered mag- 
netization does not influence superconductivity, however 
large it is. 

The heavy fermion superconductors, such as UPd2Al3 
and CePt 3 Si, can be candidates of the present mecha- 
nism. The antiferromagnetic long-range order is con- 
sidered to be characterized by the wave vector Q = 
(0,0,tt/c), both in UPd 2 Al 3 Qj and in CePt 3 Si 0. Su- 
perconducting transitions have been observed at T c = 
2.0 K and 0.7 K, below the antiferromagnetic transition 
temperatures Taf = 14.3 K and 2.2 K, in UPd 2 Al 3 
and CePt 3 Si Q, respectively. It has been suggested 
that the magnetic moment is large, i.e., 0.85/is/U.in 
UPd 2 Al 3 0, but small, i.e., 0.16^ B /Ce, in CePt 3 Si 0. 

The order parameter of interlayer spin-singlet pairing 
has a horizontal line node. This also agrees with the 
experimental results in the compound UPd2Al 3 . The ex- 
istence of the line node is suggested by the nuclear mag- 
netic resonance (NMR) measurement [5|, in which the 
Hebel-Slichter peak was absent, and T\~ oc T 3 was ob- 
served. The singlet state is supported by the presence 
of NMR Knight shift y|, and the Pauli limited upper 
critical field . In recent angle resolved magnetother- 
mal transport measurements 8], two- fold oscillation in 
the rotation perpendicular to the ab plane was observed, 
while no oscillation was observed in rotation in the ab 
plane. These experimental results are reproduced by the 



order parameter of the form of A(fe) = Ao cos(fcjc). 

Coexistence of superconductivity and magnetism has 
been studied in various models by many authors 0, 
ITlj . In particular, spin singlet superconductivity in the 
presence of ferromagnetic layers has been studied by 
many authors 0. In the models examined in those pa- 
pers, superconductivity occurs in a subsystem different 
from the magnetic layers. In contrast, we examine a 
model in which superconductivity occurs in electrons on 
magnetic layers in the present model. Kopaev also stud- 
ied superconductivity when only magnetic electrons are 
present [lOj , although it exists only in the vicinity of the 
domain wall. 

The magnetic structure mentioned above can be mod- 
eled most simply by the Hamiltonian 



H — Hq + Hjj + Hj 
with the kinetic energy term 

Ho = ^2^ k ) c ia c ^, 

ha 

the on-site Coulomb interactions 

i 

and the exchange interactions 



(1) 



(2) 



(3) 



(4) 



c}~Cia, where <x denotes the vector of Pauli 



We have defined Si — h V . ct 
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and m a = c\ a 
matrices, and Cka and Ci a denote the electron operators. 
We define J tj = J > for R 3 = R t ± c, J tj = - J\\ < 
for nearest neighbor sites on the same layer, and 
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otherwise. Here, Ri and c denote the lattice 



vector of the site i and the unit lattice vector in the c- 
direction. The interlayer antiferromagnetic exchange in- 
teraction J originates from the interlayer superexchangc 
or kinetic exchange processes of electrons. The intralayer 
ferromagnetic exchange interaction Jm expresses the ef- 
fect of the exchange Coulomb interaction, which is usu- 
ally smaller than U and J, but necessary to stabilize the 
present magnetic structure. The ferromagnetic correla- 
tion in each layer is due to U and J|| . However, two- 
dimensional (2D) long-range order without order in the 
c-dircction cannot occur due to the thermal fluctuations 
in the present isotropic model, however large U is. Tran- 
sition to the long-range order occurs only in the presence 
of interlayer exchange interaction J, and the transition to 
the three dimensional antiferromagnetic long-range order 
at T = Taf is the only magnetic transition. 

Many examples of compounds which can be modeled 
by 2D Heisenberg ferromagnetic model with antiferro- 
magnetic interlayer exchange interactions are summa- 
rized by a review article by Jough and Miedema |l2j | . For 
example, it was obtained from experimental data that 
J/J,, w 8 x 1(T 3 , 3.4 x 1CT 3 and 0.21, Taf = 13.8, 
16.8 and 18 K, and J\\/k B = 18.8, 5.25 and 3.0 K in 
the compounds Rh^CuCLi, CrCi3 and NaCrS2, respec- 
tively. These compounds have the ferromagnetic short- 
range order in each layer at temperatures higher than 
the transition temperature due to the intralayer ferro- 
magnetic exchange interaction Jy, and undergo the long- 
range order by the weak interlayer exchange interaction 
J at the transition temperature Taf- In many exam- 
ples, the interlayer exchange interactions are antiferro- 
magnetic and much weaker than the intralayer interac- 
tion in most cases. In our model, we also take into ac- 
count the on-site Coulomb repulsion U in addition to the 
intralayer exchange interaction J|| to stabilize the fer- 
romagnetic structure in each layer. Later, we consider 
a situation in which U 3> Jii as an example, but it is 
straightforward to apply the theory to the opposite case 
U < J,,. 

It is well-known due to the Mermin- Wagner's theo- 
rem [Hi ] that purely 2D isotropic Heisenberg model can- 
not exhibit any long-range order at any finite tempera- 
ture. For the long-range order to be stabilized, an ad- 
ditional Ising type intralayer interaction or a three di- 
mensional (interlayer) interaction is necessary. However, 
the former does not stabilize the present antiferromag- 
netic configuration in the stacking direction, as observed 
in UPd 2 Al 3 and CePt 3 Si 0. Therefore, more or 
less, interlayer antiferromagnetic interaction must exist 
in the present compounds. The physical origin of the an- 
tiferromagnetic interlayer interaction is interlayer kinetic 
exchange or superexchange process. In the former pro- 
cess, the interlayer exchange interaction J is written as 
J ~ 4t\/U , where t± denotes interlayer electron hopping 
energy. Since t± is expected to be small from the crystal 
structure, the perturbation theory to derive the above 
experssion of J would be justified. 



When we apply the Hamiltonian eq. to the com- 
pound UPd2Al 3 , we should note that it has been sug- 
gested by thermodynamic measurements that the 
magnetic and superconducting transitions occur in nearly 
disjunct subsystems in this compound. However, even if 
this is true, the present theory holds if the superconduct- 
ing subsystem has a similar interlayer exchange interac- 
tion, which is plausible because both of the two subsys- 
tems have a 5/ character and coexist in the same crystal 
structure. We discuss an application of the present the- 
ory taking into account the two-fluid model later. 

The interaction terms can be rewritten as 

Hu + Hj = — ^ v ( k > fc '> <?) 4+qt Cfc T c t-qi Cfc'l ( 5 ) 



kk' q 



with 



V(k, k' , q) = U - J(q) - J(k - k! + q) 
+ J {l (q) + J l] (k-k' + q), 



(6) 



where J(q) = J cos(q-c) and J\\{q) are the Fourier trans- 
forms of the interlayer and intralayer exchange interac- 
tions, respectively. Since we have not specified the lat- 
tice structure of the layer, the expression of J11 (q) is not 
shown, but it does not depend on q z and must have a peak 
around q\\ — 0, where q\\ = (q Xl q y ). Similarly, the form 
of £(fc) also depends on the lattice structure of the layer. 
We have simplified it as £(fc) — £,\\(k\\) + £±(fc z ) with 
£||(fc||) = H 2 \k ll \ 2 /2m* - [i and fx(M = -2t ± cos(k z c) 
for convenience, where fc|| = (k Xl k y ). This simplification 
does not essentially change the qualitative results. In 
this paper, we examine the system in which t± and Jy 
are small. We take units with h = ks = 1. 

First, we describe the magnetic transition. Let us ex- 
amine the spin propagator 



X(Ri - Rj,r) 



-<t t s?(t)s;(o)) 



(7) 



in the random phase approximation (RPA). We define 
the Fourier transform %(q, iv m ) by 



X(q,i"n 



i(q-R l —i' rn r) 



x(Ri,r), (8) 



where v m = 2irmT denotes the Matsubara frequency, 
and j3 = 1/T. If x(Q,®) diverges, it indicates the phase 
transition to the magnetic long-range order with Q. If 
we omit t± and Jm , we obtain 



1 



Xn(q-iv-m) 



1\-\U- J(q)}xo(q,iv m )' 



(9) 



where 



k n (10) 
x G { °\k + q,iuj n + iv m ) 
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with the bare electron Green's function G^ ' . When t 
0, the free susceptibility xa is expressed as 



Xo(g,0) = po 1 - Re 




(11) 



at T = 0, where po = m*ab/2nh 2 denotes the density 
of states of the 2D system. We have defined the effec- 
tive mass m* , the in-plane Fermi momentum kp\\, and 
the lattice constants a and b. The maximum xi^i^m) 
occurs at arbitrary q — (q\\,ir/c) with |qiy| < 2fejrii and 
v m = 0. This degeneracy is removed by Jii ^ 0, which 
is small but exists in practice. Hence, xilt^m) reaches 
its maximum at q = (0, 0, tt/c) = Q. It is easily verified 
by replacing U — J with U — J + Jii in eq. © . Further- 
more, when we take into account t± ^ 0, xo has a peak 
around q = (qrii , 7r/c) = Q' with |gii 
difference Xo(Q',0) - Xo(0,0) oc t\ 
when tj_ ^ 0, we must assume Ji 

but sufficiently large for the maximum of x to occur at 
q = Q, so that the magnetic order of Q is stabilized. 

When these conditions are satisfied, antiferromagnetic 
transition occurs at a temperature which satisfies 



2kpu, but the 
is small. Therefore, 
I 7^ which is small 



! = (£/ + J) xo(Q,0) 



(12) 



from eq. ©, where X o(Q, 0) = Xo(0, 0) = po/{er^ + 1). 
The chemical potential \i is determined by the equation 



2/?-Voln(l + e^). 
-0n/2 Pa y Therefore, 

we obtain the antiferromagnetic transition temperature 



for the electron number per site n 
Thus, we obtain xo(Q, 0) = po(l — 



Taf = 



(U + J)po 
(U+J)p -l 



(13) 



2po In 



to an ordered state with the wave vector Q = (0, 0, tt/c) 
when (U + J)po > 1, while Taf — otherwise. 

In the antiferromagnetic phase, the electron states are 
affected by spontaneous staggered magnetization. We de- 
fine A and B sublattices (sublayers) whose majority spins 
are up and down, respectively. We write the electron op- 
and bj a for i 6 4 and j G B. Therefore, 



erators as di C 
we have 



{a\ a a ia ) 



(14) 



am, 



where {Sf) = —{Sj) = m for i G A and j G 5. We have 
defined <r = +1 and —1 in equations, which correspond 
to a =f and J, in suffixes, respectively. Corrections to the 
kinetic energy due to Hjj + Hj are taken into account by 
the mean field approximation as 



Hmf = 



ka 



OtlMF 



[a{ a a ka - bl a b ka ] 



(15) 



with Hmf — (U + J + z» J\\/2)m, where z\\ denotes the 
number of nearest neighbor sites in the layer, and the 



summation Y^'k 1S carried out over the half Brillouin zone. 
Therefore, the total kinetic energy term Ho = Hq + Hmf 
is written as 



U(al a b ka 



(16) 



where = ^ —gIimf and £^ = £|| +<tIimf- The mean 
field approximation is consistent with the RPA, which we 
have used to derive Taf, although we have neglected Jy 
in eq. l|T3)l. 

Now, let us examine superconductivity. We will show 
its formulation in the antiferromagnetic phase, but it is 
immediately reduced to that in the paramagnetic phase 
by putting m = 0. The exchange interaction Hj con- 
tributes to pairing interaction, while it causes the anti- 
ferromagnetic transition and creates the exchange field. 
We rewrite Hj as 



Hj 



(17) 



where = 2" 

■M 



1 / 2 (ai|6 J | — aiibji). The statistical av- 
erage i^fj'} is the order parameter of interlayer spin- 
singlet pairing. Here, we have neglected Jy, since it does 
not have an important effect on superconductivity if it is 
small. In eq. Q17|). it is found that J contributes only to 
spin-singlet pairing as an attractive interaction in its first 
order. In the BCS approximation, eq. 1)17(1 is written as 



Hj W ^ 0"[A(fe) a k ab-k-a + tl.C.] 



with the order parameter 

A ( fe ) = -^fE'^ - k ') < fo V CT <«L<>- 



2N- 

k' a' 

Therefore, we obtain 



A(fe) = Aq cos(fc z c) 



with 



J 



A ° = -^?E v'cos(k' z c)(bl k ,_ a ,a 



2N^ 

k'cr 



h' a 1 



(18) 



(19) 



(20) 



(21) 



In the same approximation, the on-site Coulomb interac- 
tion H\j is ineffective for anisotropic superconductivity. 

It has been proposed that in the higher order of J, 
the pairing interaction is enhanced by the exchange of 
magnons |16| . This mechanism has also been examined in 
the compound UPd 2 Al 3 [H El El HI 

However, since 

the spin fluctuations are weak at temperatures much 
lower than the antiferromagnetic transition temperature 
(T < Taf/7), the pairing interaction mediated by the 
magnons is weak [l6j. Hence, in this paper, we neglect 
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them in comparison to the direct pairing interaction de- 
scribed in eqs. (|T7|l and lfT5)l. The present direct pairing 
interaction is not mediated by the magnons. In a broader 
sense, however, one may regard the present pairing inter- 
action as mediated by the spin fluctuations, because the 
superexchange interactions are derived from the virtual 
process of electrons, with which their spins correlate. 

When J = 0, the present model is reduced to the quasi- 
2D Hubbard model. In the perturbation theory based on 
it, it was shown that antifcrromagnctic fluctuations in- 
duce intralayer singlet pairing near the antifcrromagnetic 
phase [2l|. However, the present system is ferromagnetic 
in each layer. In the absence of J, the propagator of 
the fluctuations x{q) has a broad peak around q« w 
as we can see in eq. (0, in contrast to the sharp peak 
in the antiferromagnetic case, where the Fermi-surface 
nesting occurs. Therefore, the pairing interaction is not 
strongly enhanced by the spin fluctuation unless U is 
large. Nisikawa and Yamada [2^1 examined the UPd2Al3 
on the basis of a 2D Hubbard model taking into account 
the lattice structure, although they did not examine in- 
terlayer pairing. In the presence of interlayer interac- 
tions, t± and J, the spin fluctuations will enhance the 
interlayer singlet pairing interaction. 

From eqs. (|16H and (|18fl . we obtain 



(22) 



where we have defined ip^ = (aj^ a,-ki fr-fcx) and 

/ ef T a o -A* \ 

M(k)= IIT 

-A -U 

The Green's function is defined in matrix form by 
From the equation of motion, we obtain 

G(k,iu n ) = [iloJ -M{k)Y\ 



(23) 



(24) 
(25) 



where / denotes the 2x2 unit matrix. We obtain 
the quasi-particle energies, ±Eji.(k) and ±Es(k), where 

7'2 _ £2 _l < 'I : 1,2 

|A| 2 ^2/ 2 ]i/2. if W e define G t 



E% = E 2 + F, E%=E 2 - F, E*=a+e ± + h\ 1F + | A| 



and F = 2^1 + ^) 
as the component of the matrix G, we have 

^tT a - fe 'l)= G 32(fc , ,T = -0), 

(6 t _ fe , i 4, T ) = -G 4 i(fe , ,r = -0). 
Therefore, we obtain the gap equation 

1 i tanh Ex{k "> 

A(*)-^(*-«0E 



(26) 



X=A,B 



-22_A(fc'), (27) 

2E x (k') V h y ' 



FSlA 



FSTA 




FS i B FS T B 



FIG. 1: Schematic figure of the Fermi surfaces of J and J. 
electrons in the A and B layers, and interlayer spin-singlet 
pairing. The abbreviation FSaX denotes the Fermi surface of 
the spin a electrons on the X sublattice, where a =|, J, and 
X = A, B. The splits of the Fermi surfaces of the | and J. 
spin electrons do not affect interlayer spin-singlet pairing. 



In the limit of Ao — > 0, we have Ex 

f||+£x, £,b = £|| -Sx, and <5_l = \/h 
we obtain the equation for T c 



2 

MF 



|gx|, where £ A = 
+ Therefore, 



tanh[£x(fc')/2T o ; 



AN 



X=A,B 



(28) 



It is easily verified that the right hand side of eq. 128|) 
exhibits a logarithmic divergence in the limit of T — > 0, 
irrespective of the value of the exchange field /imf- Loga- 
rithmic divergence, which results from the Fermi-surface 
instability, is crucial for the occurrence of superconduc- 
tivity. Since the divergence occurs even if Iimf is very 
large, the existence of a large local magnetic moment as 
observed in UPd2Al 3 does not deny the occurrence of 
superconductivity in the electrons which have the present 
magnetic structure. It is interesting that superconduc- 
tivity is not influenced by large spin polarization in each 
layer, although such polarization creates the strong ex- 
change field on the electrons responsible to superconduc- 
tivity. It is well-known that in ferromagnets, strong ex- 
change field suppresses superconductivity by Pauli para- 
magnetic pair-breaking effect. 

This result can easily be verified as follows. The 
present superconductivity is due to singlet pairing of aka 
and b-k-a electrons, i.e., spin a electrons on A sublattice 
and spin —a electrons on B sublattice. When we define 
the Fermi surface of each sublattice, the magnitudes of 
the Fermi momenta of those electrons are equal, irrespec- 
tive of the magnitude of the Zeeman splitting due to the 
exchange field in each layer, as schematically shown in 
Fig. ^ Therefore, the present pair states are not influ- 
enced by the magnetic moments. 

Needless to say, interlayer pairing does not mean that 
the coherence length in the c-direction £oi- is on the order 
of the layer spacing. In the present system, £o_l is on 



5 



the order of Vp j_/Aq, where Vp± and Ao denote Fermi 
velocity in the c-direction and the scale of magnitude of 
the order parameter at T = 0, respectively. We obtain 
£o± > c, if t± > A . 

When tj_ is negligible, the transition temperature T c is 
obtained as follows. Since S± = |/imf|, we can integrate 
cos 2 {k' z c) with respect to k' z first in eq. I|28p. Unless S± > 
p, the density of states is constant when t± « 0. The 
fcy integral is approximated by the £y ± Hmf integrals 
with an effective cutoff energy W c , which is on the order 
of the band width. More explicitly, it is expressed as 
W c = \{W - /x - h M F)(W + h M F){p + h M F)(v - 
/imf)] 1 ' 4 , where W and p denote the band width and 
the chemical potential measured from the bottom of the 
band, respectively. Carrying out the integral, we obtain 

T c = 1.13W c e~ 2/Jpa . (29) 

Here, it is found that T c is not influenced by the sponta- 
neous staggered magnetization j/ism, as expected from 
the above argument. 

Figure |21 depicts the antifcrromagnetic transiton tem- 
perature Taf and the superconducting transition tem- 
perature T c scaled by n/p and W c , respectively. These 
scales have the same orders of magnitude, i.e., po ~ 
1/W C . It is found that, as the order of the magnitudes, 
the experimental results in UPd2Al3, Taf — 14.3 K and 
T c = 2.0 K could be explained within the present mech- 
anism, if the effective band width is on the order of 
W = 10 ~ 100 K, which is realistic for heavy fermion 
systems. 

Figure |21 shows the values of Jpo and Upo for a given 
ratio of a = (T AF po/n)/{T c /W c ) = T AF /T c xW c p /n. If 
we consider W c po/n ~ 1, the experimental data Taf = 
14.3 K and T c = 2.0 K give a ~ 7. Therefore, the exper- 
imental value of ratio a can be reasonably reproduced 
for moderate values of the coupling constants J and U. 
We will compare the theoretical and experimental results 
more closely below. We obtain Taf > T c or Taf = 0, 
except in a very small region of the phase diagram. 

We note that the resultant singlet order parameter is 
invariant under rotation in spin space. The rotational 
transformation is made by 

6 6 9 

R(9) = exp[i— <jy\ — cos — + ia y sin — . (30) 

The electron operators in the rotated space are defined 
by 

it) Mi} (3i) 

where Ck a — Ofco- or b^a- It is easily verified that 

T,^ b - k -A,) = J2^-A a )- (32) 

cr a 

In addition, the present pairing interaction is rotation- 
ally invariant, i.e., J(k,k') does not have spin suffixes. 



1| . . . . , . . . . n1 




FIG. 2: The transition temperatures as functions of Jpo- 
The solid curves (a) - (d) show the results of poTaf /n for 
Upo = 0.2, 0.5, 1.0, and 1.5, respectively. The dashed curve 
shows the results of T c /W c - 
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0.5 1 



FIG. 3: Contour lines for various ratios Taf/T c in J-U 
plane. The solid curves (a) - (f) are the countour lines 
for a = 2, 3, 5, 10, 20, and 100, respectively, where a = 
{Taf po/n)/ {Tc/Wc) = T AF /T C x W c p /n. The dotted curve 
is that of a = 7, which corresponds to the experimental result 
of UPd2Al3. The dashed curve shows the phase boundary be- 
tween the antiferromagnetic and paramagnetic (PM) phases 
at T = 0. 

Therefore, A(fc) given by eq. i|19|) is isotropic in spin 
space, irrespective of the direction of the magnetic or- 
der. In fact, the resultant gap equations l|27|) and l|28|) 
are invariant under spin rotation. 

These results may explain the experimental data of 
the muon spin rotation measurements in UPd2Ai3 |23| . 
in which it was observed that the London penetration 
depth and the magnetic susceptibility reduction below 



6 



T c are essentially isotropic, while the total susceptibility 
remains strongly anisotropic. In the present mechanism, 
not only the singlet nature but also the rotational in- 
variance of the pairing interaction play essential roles in 
the spin isotropy of A(fc). In contrast, in the "magnetic 
exciton" mechanism, the pairing interaction and the or- 
der parameter are anisotropic if the magnetic system is 
anisotropic. 

The compound CePtaSi is another candidate for the 
present pairing mechanism, although the ratio of the 
transition temperatures Taf/T c is much smaller than 
that in UPd2Al3. In contrast to UPd2Al3, the critical 
field largely exceeds the Pauli paramagnetic limit Hp 
estimated by the simplified formula Hp « 1.86[T/K] x 
T C [K] w 1.4T in CePt 3 Si @|. If the present mechanism 
of singlet pairing is realized in CePtaSi, the large critical 
field cannot be attributed to equal spin pairing. It can be 
explained by an effect of exchange field created by coex- 
isting antiferromagnetic long-range order, which reduces 
the Pauli paramagnetic pair-breaking effect |24| . It is still 
controversial whether dominant pairing in CePtaSi is of 
singlet or triplet. Recent NMR data suggests that the 
gap function may have some novel structure p5|. Ther- 
mal transport measurements have suggested that the or- 
der parameter has line nodes (2||, which is consistent 
with the present theory. It is known that the compound 
CePtaSi does not have inversion symmetry. The Rashba 
interaction j2?J has been examined to include it. Yip 
predicted that the Knight shift vanishes in the supercon- 
ductors with strong Rashba interaction. |2g 

For more close comparison with the experimental data 
of UPd2Ai3, we consider a two-fluid model 
order to take into account the model, we assume two 
different renormalization factors Z s and Z m for elec- 
trons responsible to superconductivity and antiferromag- 
netic order, respectively. We should note that in ac- 
tuality there are strong renormalization effects in the 
heavy fermion system, not only in the density of states, 
but also in the vertex corrections. It is easily verified 
by diagramatical consideration that in terms of Z s and 
Z m , the electron mass too, the density of states po, the 
band width W, and the coupling constants J and U 
are modified as rh a = Z a m , p a = Z a p , W a = W/Z a 
J a = J/Z^ and U a = J/Z%, respectively, where a = 
s,m. The experimental value 71 = 115mJ/K 2 mol 
gives p s = 1/476 K^ 1 . Equation (|29[1 is rewritten as 

T c = l.UWce" 1 ^ with A = J s p s , p s ~ 1/W S , and 
W c = W s /2. Hence, T c = 2.2 K gives J s = 173 K. The 
antiferromagnetic transition temperature of eq. (|13fl is 
written as 



~ ~ 1 ( U m -f- J ra ) pjji 
2p m In — = — 

( ~\~ Jrn)Prn 1 

Here, we simply put n — 2 for order estimations. From 
the values of J s and p s estimated above, the values of J m 
and p m can be obtained, if the ratio Z m jZ s is known. 
Therefore, we only need the values of Z m jZ a and U m for 



estimation of Tap. However, since they are not known for 
UPd 2 Al 3 at the present, we need to assume them. In the 
assumption, we require that W m is smaller than Taf con- 
sistently with the observation of a large local magnetic 
moment in UPd2Al3 |l4|. Physically, it is also plausible 
that U m is not much larger than Taf but larger than 
W m . As an example, let us assume that Z m /Z s = 30, 
which gives W m ss 8K < T A f and J m « 0.19 K. In 
this case, if we assume U m — 20 K and 30 K as exam- 
ples, we obtain Taf ~ 10 K and 21 K, respectively. As 
another example, we assume Z m /Z s = 20, which gives 
W m w 12 K < T AF and J m w 0.43 K. In this case 
U m = 30 K and 40 K give T AF ~ 16 K and 27 K, respec- 
tively. The values J m w 0.19 K and 0.43 K obtained in 
these examples are not outrageous as energy parameters 
in real materials. In fact, J fa 0.15 K and 0.63 K were ob- 
tained from experimental data in Rb2CuCl4 and NaCrS2, 
respectively [l2|- Here, J was not estimated as the bare 
parameter but estimated as the dressed (observed) pa- 
rameter like J m . Although our estimations are crude, 
the values of Taf obtained above are on the same or- 
der of the experimental result Taf = 14.3 K. Therefore, 
we find that the present mechanism reproduces consis- 
tent orders of magnitudes of T c and Taf for appropriate 
values of Z m jZ s and U m . 

The result that Taf > T c except for in a very limited 
region unless Taf can be explained physically as fol- 
lows. In the present mechanism, both interlayer antifer- 
romagnetic long-range order and interlayer singlet super- 
conductivity are induced by interlayer antiferromagnetic 
exchange interaction. However, strong on-site repulsion 
contributes to stabilization of the ferromagnetic structure 
in each layer, while it does not contribute to interlayer 
singlet pairing. Therefore, the magnetic transition oc- 
curs at a higher temperature than the superconducting 
transition temperature unless the on-site U is negligibly 
small. 

Interlayer pairing has been studied by many au- 
thors 29]. In this paper, we have examined the mag- 
netic mechanism of the pairing interactions for interlayer 
singlet pairing, when the electrons are on the magnetic 
layers. However, irrespective of the pairing mechanism, 
interlayer pairing of the present type seems to be the 
most favorable, apart from equal spin pairing, when the 
present type of antiferromagnetic long-range order co- 
exists. Other pairing states, such as intralayer singlet 
pairing, are strongly suppressed by the splitting of the 
Fermi-surfaces of the electrons with up and down spins 
due to the antiferromagnetic moment. Even in the two- 
fluid model, the exchange field must be induced on the 
electrons responsible to superconductivity. 

In conclusion, in antiferromagnets with the magnetic 
order of the wave vector Q = (0, 0, 7r/c), magnetic inter- 
actions may induce the superconductivity of interlayer 
spin-singlet pairing, the order parameter of which has a 
horizontal line node. It was found that superconductiv- 
ity and an antiferromagnetic long-range order with large 
localized magnetic moments to can coexist, and that T c 
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is not influenced by the magnitude of m. It was also 
found that Taf > T c in most cases, unless Taf = 0. 
The present model may describe an essential aspect of 
antifcrromagnetic heavy fermion superconductors, such 
as UPd2Al3 and CePtsSi. The orders of the magnitude 
of Taf and T c and their ratio Taf/T c « 3 ~ 7 can be re- 
produced by assuming moderate parameter values. The 



resultant order parameter is consistent with the observa- 
tions mentioned above 0, H SH IH El • 
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